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Abstract 

Using the symplectic tomography map, both for the probability 
distributions in classical phase space and for the Wigner functions of 
its quantum counterpart, we discuss a notion of Lyapunov exponent 
for quantum dynamics. Because the marginal distributions, obtained 
by the tomography map, are always well defined probabilities, the 
correspondence between classical and quantum notions is very clear. 
Then we also obtain the corresponding expressions in Hilbert space. 

Some examples are worked out. Classical and quantum exponents 
are seen to coincide for local and non-local time-dependent quadratic 
potentials. For non-quadratic potentials classical and quantum expo- 
nents are different and some insight is obtained on the taming effect 
of quantum mechanics on classical chaos. A detailed analysis is made 
for the standard map. 

Providing an unambiguous extension of the notion of Lyapunov 
exponent to quantum mechnics, the method that is developed is also 
computationally efficient in obtaining analytical results for the Lya- 
punov exponent, both classical and quantum. 

*on leave from the P. N. Lebedev Physical Institute, Moscow, Russia 
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1 Introduction 



Classical chaotic motion is characterized by the existence of positive Lya- 
punov exponents or positive Kolmogoroff-Sinai entropy. Because the defini- 
tion of these quantities is based on the properties of classical trajectories in 
phase space, it is not obvious what the corresponding quantities in quantum 
mechanics should be. This situation led to the proposal of many differ- 
ent quantities to characterize the quantum behavior of classically chaotic 
systems (|TJ - To be sure that one constructs quantum mechanical 

functionals, with exactly the same physical meaning as the classical quan- 
tities which characterize classical chaos, the natural suggestion would be to 
use also a phase-space formulation for quantum mechanics, rather than the 
usual Hilbert space formulation. The difficulty here lies in the fact that 
quantum phase-space is a non-commutative manifold with the usual point- 
wise product of functions being replaced by the *-product[^7}]. One possible 
solution is to use the tools of non-commutative geometry for this formulation. 
Another approach, however, is to look for (commutative) quantities which 
have the same formal structure both in classical and in quantum mechanics. 



The Wigner function p8[, which some authors have attempted to use for this 
purpose, is not the appropriate choice because, unlike the classical proba- 
bility distributions, it is not positive definite. In fact, the Wigner function 
is only correctly interpreted in a non-commutative geometry setting or, al- 
ternatively, as an operator symbol which by the Weyl map corresponds to a 



well defined operator in Hilbert space P7 



There is however a set of phase-space quantities that have the same math- 
ematical nature in both classical and quantum mechanics. This is the set of 
marginal distributions of the symplectic tomography formulation ( ||29|| - ||36||), 
which are in both cases well-defined probability distributions. In Section 2, 
we review the symplectic tomography formulation of classical and quantum 
mechanics. In both cases the dynamics is defined by a set of marginal proba- 
bility distributions. The difference between classical and quantum mechanics 
comes only on the modification of the equations of motion, which in this for- 
mulation is the analog of the Moyal def ormat ion |37|| p7 of the phase-space 
algebra. 

Once the appropriate phase-space quantities are identified and classical 
Lyapunov exponents (and local entropies) are formulated in terms of proba- 
bility distributions, the transition to quantum mechanics is rather straight- 
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forward.. This is discussed in Section 3. Having obtained the Lyapunov ex- 
ponent for quantum mechanics from the marginal probability distributions, 
one is then also able to obtain the corresponding Hilbert space expression. 

Some examples are then computed, involving kicked systems on the line, 
on the 2-torus and on the circle. Classical and quantum exponents are seen 
to coincide for time- dependent local and nonlocal quadratic potentials. For 
non-quadratic potentials classical and quantum exponents are different. A 
characterization is obtained for the origin of the taming effect of quantum 
mechanics on classical chaos in the standard map. This is a step towards 
a rigorous characterization of this effect because it refers, as in the classical 
case, to the behavior of the Lyapunov exponents and not to indirect quantities 
like the energy growth. 

2 Symplectic tomography of classical and quan- 
tum states 

2.1 Classical mechanics 

States in classical statistical mechanics are described by a function p(q, p), 

the probability distribution function in 2n-dimensional phase space (q G R n ,p G R n ), 

with properties 



P{q) and P(p) being the probability distributions for position and momentum 
(the marginals of p). 

The function p (q, p) is normalized 



We consider an observable X (q, p), that is, a function on the phase space 
of the system. The inverse Fourier transform of the characteristic function 
(e tk ' x ) for any vector observable X (q, p) 







(1) 
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is a real nonnegative function which is normalized, since 

w(Y) = J p (q, p) 5 n (X(q, p) - Y) d n qd n p (3) 



and 



J w (Y) d n Y — J p(q, p) d n qd"p — 1 (4) 

As a classical analog of the quantum symplectic tomography observable, 
introduced in [p5f| , we consider the following classical observable p3[ 

X (q, p) = H ® q + v ®p (5) 

where © denotes the componentwise product of vectors 

(// © q)i = mi 

and n and v are vector-valued real parameters. Together with 
P(q, p) = - ® V®q+ (- + 1 J ®P 

v \fl ) 

(|5]) is a symplectic transformation?] of the position and momentum observ- 
ables. (1 stands for the unit vector in R n ) 

The vector variable X (q, p) may be interpreted as a coordinate of the 
system, when measured in a rotated and scaled reference frame in the classical 
phase space. For the coordinate (||) in the transformed reference frame, we 
obtain from Eq.(|2|) the distribution function (the tomography map) 

w (X, /i, v) = — ^ / e -ik'(x-„® q ~^®P) p ( ?> p ) d n qd n pd n k (6) 
[2 it) J 

This function is homogeneous 

w (XX, A/i, Xu) = \X\- n w (X, fi, v) (7) 

^^This is not, of course, the most general symplectic transformation. In general 
X l (x,p) = n\x k + v % k -p k with the corresponding expression for P(x,p). Here we have 
considered the particular case where the tensors /x and v are diagonal. That is the reason 
for the non-covariant look of our equations. 
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and Eq.(||) has an inverse 

P (?) P) — 7 — -^yn / w {X> P-, v ) ex P [* (X — p © q — v © p) • 1] d n X d n p, d n v . 
(An 2 ) J 

(8) 

Since the map 

P (q, P) =>w A*, v) 

is invertible, the information contained in the distribution function p (q, p) is 
equivalent to the information contained in the marginal distributions w (X, fi, v). 

The Boltzman evolution equation for the classical distribution function 
for a particle with mass m = 1 and potential V(q), 



d^df ^ + P • %P (g, P, t) - V q V(q) • V pP (q, p, t) = 
can be rewritten in terms of the marginal distribution w (X, fi, u, t) 



(9) 



— -fi*V u w- V X V (-V^ 1 © Vj • (z/ © V x w) = (10) 
For the mean value of the position and momentum we have 
q\\ f -i- -\ ( 1 



pJ/ f P(Q,P)[1 j d n qd n p (11) 

= i Jw(X,fi, u)e iXml ( ^ J (6 n (p,)5 n {v)) d n Xd n pd n v 

With a change of variables X — > p,X and the homogeneity property w(fiX, //, 0) = 
fi~ n w (X, 1, 0) one may, for example, check the consistency of this definition 
for the mean value of the position 

(q) = J w(X,l,0)X<TX (12) 

By Il cl (X, /i, u, X', /i', z/', t 2 , ti) we denote the classical propagator that 
connects two marginal distributions at different times to an d t (t > to) 

w (X, /I, u,t) = J n cl (X, //, v, X', //, v\ t, t ) w (X', /i', v', t ) d n X' d n p! d n v> '. 

(13) 
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The propagator satisfies the equation 



dIL 



cl 



- fi • v,n cl - v x v (-v^ 1 © v„) • {v © v x n cl ) = o (w) 



dt 2 

with boundary condition 

lim n cl (X, v, X\ n', z/, t 2 , tx) = 5 n (X - X') 5 n (// - //) 5 n {y - v') (15) 

2.2 Quantum mechanics 

For quantum mechanics the construction is similar and the mathematical 
nature of the quantities that are constructed is the same, because it is a 
general fact that the inverse Fourier transform of a characteristic function 
is a positive distribution. The marginal distributions that are obtained are 
simply related to other well known quantum mechanical quantities. It was 
shown 1 29] that for the generic linear combination 

X = fi®q + u®p (16) 

where q and p are the position and the momentum, the marginal distribution 
w (X, fi, v) (normalized in the variable X and depending on two vector- 
valued real parameters fi and v) is related to the Wigner function W(q, p). 
For n degrees of freedom one has 

/d n k d n q d n r> 
exp [-ik •(X-fi®q-u®p)] W(q, p) J F (17) 

We see that Eq.(|I?D is formally identical to (|6|) of the classical case. For a 
pure state with wave function \l/ (y), the marginal distribution would bc|35| 



w (X, /I, v) 



(2tt)" \v x ■■■v n 



r n / 2 v \ 



2 



(18) 

Eq. ([l7|) may be inverted and the Wigner function expressed in terms of the 
marginal distribution, like in the classical case of Eq. 

W(q, p) = J w i X > v) exp[i (X - /J ® q - v ® p) • 1] d n fi d n v d n X 

(19) 
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Therefore the usual quantum mechanical quantities can be reconstructed 
from the marginal distributions. These quantities (wave function and Wigner 
function) have a nature quite different from the classical quantities, however 
the marginals w (X, li, v) are in both cases positive distributions with the 
same physical meaning. 



As was shown in pi| , for a system with Hamiltonian 

2 

H = V - + V(q) , (20) 
the marginal distribution satisfies the quantum time-evolution equation 

% - Li . V v w - z\{V (-V x x © V, - i\hv © Vx) (2 , s 

-V (-V^ 1 ®V fl + i±hv®V x )}w = K ' 

which provides a dynamical characterization of quantum dynamics, alterna- 
tive to the Schrodinger equation. 



The evolution equation fl2T|) can also be written in the form 



&-H*V v w- V X V (q)m(v® V x w) 

+1 E^-ir 1 (I ) 2n+1 (y ® v -)u • • • (y ® v -)- +1 w = 



( 

where q stands for the operator 



(22) 



q = -V x ! © V 



i-i 



and a sum over repeated indices is implied. 

In Moyal's [p^1 formulation of quantum mechanics in phase-space, the 



transition from the classical to the quantum structure is a deformation of 



the Poisson algebra |27| with deformation parameter h. In the symplectic 
tomography formulation, that we are describing, classical and quantum me- 
chanics are described by the same set of positive probability distributions 
w (X, li, v), the ft— deformation appearing only in the time-evolution equa- 
tion (E2). 



For the propagator 



iv (X, Li,u,t) = J II {X, li, v, X', li', v', t, t ) w {X', li', v', t ) d n X' d n Li' d n v'. 

(23) 
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§ - n . v ji - v x v {q)*{v® v x n) 



the equation is 

an 

at 

+1 E^(-i)" +1 (f ) 2n+1 Vil (^r @ © Vx) n •■■(,© v*) Wl n 

= 

(24) 

with boundary condition 

lim IT (X, /I, u, X', fi', u', t, t ) = 5 n {X - X') 5 n {p - p') 5 n {y - v') (25) 

3 Lyapunov exponents 

3.1 Density formulation in classical mechanics 

Lyapunov exponents and other ergodic invariants in the classical theory are 
usually formulated in terms of quantities related to trajectories in phase- 



space, like tangent maps, refinement of partitions, etc.|38|]. Here, as a prepa- 
ration for the formulation of Lyapunov exponents in quantum mechanics, 
using the marginal distributions w (X, fi, v) , we explain briefly how these 
quantities may, in classical mechanics, be expressed as functionals of phase- 
space densities rather than in terms of trajectories. For more details we refer 
to p. 

A density in phase-space is a non-negative, normalized, integrable func- 
tion, the space of densities being denoted by D 

D = {peL 1 :p>0,\\p\\ 1 = l} (26) 

D is the space of functions that, by the Radon-Nikodym theorem, character- 
ize the measures that are absolutely continuous with respect to the underlying 
measure in phase-space. However, to define Lyapunov exponents by densi- 
ties, it is necessary to restrict oneself to a subspace of admissible densities 
defined as follows: 

To each p e D we associate a square root, that becomes an element of an 
L 2 space. We then construct a Gelfand triplet 

E* D L 2 D E (27) 

where E is the space of functions of rapid decrease topologized by the family 
of semi-norms Ha^c^/I^ and E* is its dual. Because E is an algebra / G E 
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implies f 2 G E. Therefore for each / such that 



1 , P = P 



is an 



admissible density. The restriction to such a subspace of admissible densities 
is necessary to be able to define Gateaux derivatives along generalized func- 
tions with point support. Gateaux derivatives along derivatives of the delta 
function play for densities the same role as the tangent map for trajectories. 



In this setting the Lyapunov exponent is[12 



A„ = lim - log 

t— >oo t 



v l D aA / dp{y)yP t p{y) 



(28) 



v G R 2n , ||-|| is the vector norm and the Gateaux derivative Dq.s x operates 
in the argument of the functional, that is, on the initial density 



D 8iSx F (p(y)) = lim - {F (p(y) + ed^y — x)) — F (p(y))} (29) 

e— >0 £ 

p is the invariant measure in the support of which the Lyapunov exponent 
is being defined and P t is the operator of time evolution for densities 



P t p(y,0) = p(y,t) 



(30) 



A simple computation shows that the expression (|28|) is equivalent to 
the usual definition of Lyapunov exponent in terms of trajectories and the 
tangent map 

A„ = lim -1okIIZ>15;|| (31) 



lim - log \\DT*v\\ 

where DT l x stands for the evolved tangent map applied to the vector v at the 
phase-space point x. Here and in Eq.(^8|) x and y are phase space vectors, 
that is, in the notation of Sect. 2, x = (q x ,Px) and y = (q y ,p y ). 

According to Oseledec theorem J39[ fjOfl , for p— almost every point x there 
is a decreasing sequence of vector spaces 



R 



2n 



Ei(x) D E 2 {x) D ••• D E T = {0} 



such that, by choosing the vector v in E s (x) \ E s+ i(x), the sth Lyapunov 
exponent is obtained by the above calculation. 

A similar construction is possible for the metric entropy. For the en- 
tropy, the notion that seems more appropriate for generalization to quantum 



mechanics jIT]] |T2| is the Brin-Katok local entropy |38||, which for the classical 
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case and for a compact metric space is equivalent to the Kolmogoroff- Sinai 
entropy. It is denned as follows: 
Define 

B £ (T, t, x) = {y : d (T T (x), T T (y)) < e ,0 < r < t} 

B £ (T,t,x) is the ball of phase-space points around x that, in the course of 
time evolution, do not separate by a distance larger than e up to time t. 
T T (x) is the image of x after the time r and d(-, ■) is the distance. The local 
entropy h(T, x) measures the weighed (in the /i-measure) rate of shrinkage 
in time of the ball B £ (T, t, x), namely 

h(T, x) = lim lim { — log /x (B e (T, t, x)) 

As in the case of the Lyapunov exponent, this quantity may be expressed as 
a functional of (admissible) densities by rewriting the ball B £ (T, t, x) as 

B E (T,t,x) = {y: D {Sx _ Sy) (j dfi(z)zP T p(z^j <e;0<r<t| 

3.2 Classical and quantum Lyapunov exponents by marginal 
distributions 

Let us now translate the equations of the preceding subsection in the to- 
mographic framework discussed in Section 2. Initial densities are, by the 
tomographic map, mapped to initial tomographic densities by (|(|) 

p(q,p) ->w(X,fi,u,t = Q) = w(X,fi,u) (32) 

To compute the Gateaux derivatives notice that the generalized density (e 
E*) 

(vi *V q + v 2 * V p ) {5 n (q - g ) 5 n (p - Po)} (33) 
is mapped to the tomographic generalized density w v (g E*) 

w v (X, fi, v) = ((« i © fi + v 2 © v) • V x ) S n (X - fxq - up ) (34) 
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According to Eq.(^), to compute the Lyapunov exponent one has to obtain 
the expectation value of a generic phase-space vector on the time-evolved 
perturbation of the initial density fl3"Bp. Therefore 



A, 



lim - log 

t^oo t 



J d n qd n p 



V 



1C (q, p, q', p', t) Oi • S/ q ' + v 2 • V ? 



S n (q 1 — q ) S n (p' — Po)dq dp 
where fC (q,p, q',p', t) is the evolution kernel for densities 

p(q,p,t) = [ }C(q,p,q>,p',t)p(q',p') d n q> d n p> 



(35) 



Notice that in Eq . (|35"D the integration is carried over the flat phase-space 
measure d n qd n p. The result is equivalent to fl28|) for an invariant measure 
absolutely continuous with respect to d n qd n p. However the information and 
the dependence of the Lyapunov exponent on the invariant measure is carried 
by the choice of the initial point {qo,Po)- The set of Lyapunov exponents that 
is obtained by fl3"5| ) is therefore the one that corresponds to the invariant 
measure on whose support (qo,Po) lies. 

Using QTTJ) , Eq.(^) may now be rewritten using marginal distributions 



A, 



lim - log 

t^oo t 



J d n X d> d n ve iXmX 



5 n (p)5 n (u) IT cl (X, fx, v, X\ //, v\ t, 0) 



((«i ®fi' + v 2 ®i/)»Vx') 5 n (X' - fM'q - u'po) dX' n dp! n dv'' 

(36) 



where Il cl (X, /z, z/, X', //, z/', t 2 , ti) is the classical propagator defined in 

(HD-flD- 

Because (^8|) is equivalent to the usual definition of Lyapunov exponent, 
Eq. (|36|) , being equivalent to (p8|), is also a correct expression for the classical 
Lyapunov exponent. 

Now the transition to quantum mechanics is straightforward. Marginal 
distributions in classical and quantum mechanics satisfy formally identical 
expressions and have the same physical interpretation as probability densi- 
ties. The only difference lies on the time-evolution which in classical mechan- 
ics obeys Eq. flToD and in quantum mechanics the ^-deformed equation (|2l|) . 
Therefore the Lyapunov exponent in quantum mechanics will also be given 
by equation (|36|) , with however the classical propagator IT c i replaced by the 
quantum propagator II for marginal distributions, defined in (E3f) - (p5|). 
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4 Hilbert space expression for the quantum 
Lyapunov exponent 



As we will see in Sect. 5, Eq.(p6|) provides an efficient way to compute the Lya- 
punov exponent. However, for comparison with other approaches, it is useful 
to translate Eq.(|36|) in the Hilbert space quantum mechanical formalism. 

To simplify the notation, we consider n = 1, that is, a two-dimensional 
phase space. Generalization to the n— dimensional case is straightforward. To 
write the quantum Lyapunov exponent (^) in terms of the evolution operator 
acting in the Hilbert space of states, we use the following equation that 



relates the tomographic propagator to Hilbert space Green's functions 

n (X, /i, Z/, X', /i>', t) = ^j^^+y.^j G^a-y,^j 

xS (y — z — kv') exp 
dk dy dz da 



ik ( X' — X + fia — // — - — 



Then, the quantum Lyapunov exponent expressed in terms of Green's 
functions is 



A„ = lim - / da dz e- 2ipoZ G*(a,z,t) Vl 

t^oo t J 

+iv 2 {qo ~ z)G(a, 2q - z, t) 



dG 

ip G(a, 2q - z,t) - — (a, 2q - z, t) 

oz 



Using a complete set of wave functions ip n (x, t) this expression may be rewrit- 
ten 



A, 



lim-ln| [ da dz e 2ip ° z a i/j*(a, t)i[) n (z, t) 

t^oo t J — 4 

n,m 

{i[viPo + v 2 (q - z))[ijj m {a,t)ij)* m {2q Q - z,t)) 
+v 1 i> rn {a,t)i>*J l (2q - z,t)}\ 



In Refs.[0-|TT|| a quantum characteristic exponent has been defined, in 
Hilbert space, by using the time evolution of the expectation values on a 
5 —perturbed wave function. This definition may be compared with the 
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present marginal probability construction. Under a e5 (x — go) perturbation 
of the wave function, the density matrix changes, in leading order, by 

A ty(x)i/>(xf) = e {S' (x - g ) i/>*(x) + ^{x)8' {x - g ) } (37) 

On the other hand the marginal probability perturbation 

(vxfjL + v 2 y) S' (X - fiq - upo) 

induces a perturbation of the density matrix 

T* I i~t* \ / \ / ft* I rp 

ZPO I X — X II <r' I dy ~T \ _ ( ' \ r I Jy ~T dy 



Ap [x,x j = -e v I \ v x b — J + iv 2 [x - x J 5 |^g - 

(38) 

Comparing ([H]) and ( |3"8| ) one sees that, for the calculation of the Lyapunov 
exponent, they coincide on the diagonal terms (x = x ), but are different on 
the non-diagonal terms. In particular, the marginal density perturbation is 
not a pure state perturbation and cannot be reproduced by the perturbation 
of a single wave function. 

5 Examples: Kicked systems on the line and 
on the circle 

5.1 One-dimensional systems with time-dependent po- 
tentials 

We consider here one-dimensional systems with time-dependent potentials 
defined by the Hamiltonian 

H = £ + V(q,t) (39) 



For these systems, the Lyapunov exponent expression (|^) is 

(40) 



X v = lim - log 

t^oo t 



J dXdfidue iX ^ ^ j 5(ji)6(u)jF(X,n,u,t) 
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where F (X, //, u, t) is the time-evolved perturbation, namely 



F (X, n,v,t) = J n (X, //, z/, X', //, z/, t, 0) (^i/i' + v 2 z/) 5' (X' - fi'q - z/'po) dX'dfi'dv' 

(41) 

Passing to the Fourier transform 



G (k, fi,v,t) = ^- I e ikX F (X, fi, v, t) dX 

2,71 



one obtains 



lim - log 

t— »oo jf; 



lim - log 



dfxdu^ | ) )G(l,/i,z/,t) 

G( 2 )(l,0,0,t) 
^(1,0,0,^) 



(42) 



(43) 



where by and G^ 3 - 1 we denote the derivatives in the second and third 
arguments and G (k, //, v, t) is a solution of the equation 



"t"^ Z^n=l (2n+l)! r A 2 J c*.? 2 '^ 1 V ^ ife 9/j J ^ ~~ U 



(44) 



with initial condition 



ik 



G (k, n, v, t) = -— (vx/i + v 2 v) e ik ( qofM+p0 ^ 
2tx 



(45) 



Therefore, the computation of the Lyapunov exponents, both classical and 
quantum, reduces to the study of the large time limit of the solutions of 
Eq.d44]). Also the simple expression ( f43l) shows that, despite its apparently 
complex form, Eq.fl3~6|) is a computationally efficient way to obtain the Lya- 
punov exponent. 

We now study several time-dependent (kicked) potentials. 



5.2 Harmonic kicks on the line 

Here the potential is 
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2 00 

v (?) = E ( 46 ) 



7T 2 

n=— oo 



This system belongs to the class of time- dependent quadratic systems [[53 
and a solution may be found for the general case 

V (?) = a (t) t (47) 

Eq. flHl) reduces to 

9G dG . . dG 

ar-"a7 + a(f) ^ = (48) 



The ft— deformed part of Eq . fl44|) disappears and we obtain the (expected) 
result that, in this case, classical and quantum results coincide. Eq.(fJ^) has 
the solution 

G (k, v,t) = G (k, ^(e + e*) + ^(' + **) , ^ (* - **) + | (e - e*) , o) 

(49) 

the function e it) being a solution of 

" (t) +a{t)e (t) = (50) 

with initial conditions 

e(0) = l, e(0) = i (51) 

For the calculation of the Lyapunov exponent, the function G (k, -, -, 0) is the 
one given in Eq.(f45j). 



For the kicked case in (|4q) the function e (t) is obtained by establishing a 
matrix recurrence relation. Between times t n _i and t n we denote the function 
e(t) by e n (t). Then 

£n(*) = a n + b n t . (52) 
and the following matrix recurrence is obtained 



(53) 



with initial condition 

a = 1 , &o = i (54) 
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The recurrence relation (|53|) for the coefficients a n , b n yields the map 
£(Wi) ) — ( 9n+l 

£{t n+1 ) ) \ Pn+l 

for the position and momentum just after each kick 
The Floquet matrix in (|5f|) has two eigenvalues 




(55) 



Ai,: 



2ir 



4vr 2 



IT 



(56) 



Substituting in Eq.(|43|) one concludes that for z = ja/n < 4 the Lyapunov 
exponent vanishes and that for z = ja/7t > 4 there is one positive Lyapunov 
exponent 



A = In 



1 

l ~2 Z 



(57) 



This is just the classical result and also the quantum result for another def- 
inition of quantum exponent [j]] [TT|], already discussed in Section 4. The 
positive Lyapunov exponent corresponds simply to the situation where the 
Floquet operator spectrum is transient absolutely continuous [42[. However, 
the next example corresponds to a classically chaotic system which, when 
quantized, yields an absolutely continuous quasi-energy spectrum This 
suggests that it might be an example of genuine quantum chaos. The Lya- 
punov exponent analysis supports this conclusion. 



5.3 The configurational quantum cat 

The configurational quantum cat is a system with 4-dimensional phase space, 
for which the configuration space dynamics resembles the classical Arnold 
cat [53]. It describes a charged particle constrained to move in the unit square 



with periodic boundary conditions, under the influence of time-dependent 
electromagnetic pulses. It may be associated to the Hamiltonians^HJ 



H i = \v\ + \v\ + X2P1 + X1P2 y^ y 6(t- tit) (58) 



r 

n&Z 



or 

1 o 1 



Hi = -zfi + + ( x 2Pi + (xi + x 2 ) p 2 ) y*]8(t- nr) (59) 

neZ 
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A similar model may be constructed |7]] by considering only the kick part 
and defining the quantum theory directly by the Floquet operator. The deal 
with the Hamiltonians ( |58| ) or Q5§D by the tomographic formalism, we need 
to extend it to nonlocal potentials. 

Let a nonlocal Hamiltonian be written as 

H(x, p) = T(p) + V(x) + I(x, p) (60) 

where T(p) is the kinetic energy, V(x) the local potential energy and I(x,p) 
the symmetrized position and momentum- dependent interaction. Then, the 
equation for the density operator p 

p + t(Hp- pH) = (61) 
becomes, in the tomographic representation 
w(X,f2,v,t) +iif(x (1) ,p (1 V(X,/2,z7,t) -itf(x (2) ,-p (2 V(X,/I,z7,t) = 

(62) 

where, for n degrees of freedom X = {X\, X 2 , . . . , X n ), ft = (/i 1; p 2 , ■ ■ ■ , Mn)> 
V = (z/l, v 2 , ■ ■ ■ , *v), and the components of the vector-operators xW>( 2 ), 
pW'( 2 ), act on u>(X,/2, v,t) as follows 

_ ( 9 V 1 9 i_ _d_ _ ^ 2 )* 
\dX k J dp k + 2 Vk dX k Xk 

(63) 

.(i) _ _i_ _d / d V 1 g _.( 2 )* 

For the propagator II (X, fi, u, X', fi', z/, t) corresponding to Eq.(|62|) one 
obtains an equation where the quantum contributions are explicitly expressed 
by a series in powers of h. To do this, we first introduce some notation. 

Let n- vectors x and p be described by one 2n- vector Q a , a = 1,2, ... , 2n, 
with components (pi,p 2 , ■ ■ .,p n ,Xi,x 2 , ■ ■ ■ ,x n ), i.e., 

Ql— Pi) Q2 = P2 , Qn = Pn , Qn+1 = %1 , Q211 = X n 

Let also define the operator- vector Q with components Q a (a = 1, 2, ... , 2n) 
<9 \ 1 <9 ~ / 9 \ 4 9 ~ / <9 \ 1 <9 



Ql [dXj dv x ]Q2 \dX 2 ) du 2 ] --- ]Qn \dX n J dv n 

Qn+l — — 7; j Qn+2 = ~ 77T7 - 7j j • ■ ■ j Q2n 



dXj <9/V \dX 2 J dfi2 \9X n 
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and a 2n-vector d with the components 



2 9X X 2 <9X 2 ' 2 ^ 

, _ v\ d v 2 d v n d 

"n+1 — "7T ox^ ' — TTTTTT 5 • • ■ «2n 



2 9Xi ' ^ 2 <9X 2 ' ZK 2 dX„ 

Then, the equation for the propagator of Eq. (|5"2"|) is 



n - j| sin (mA| if (Q^II = , (64) 

where 

d dQ = T> da ao:- 

Using the series expansion for sin a and separating the classical-limit term 
from the quantum corrections one obtains 

• (9 2 °° (— l) n ( (9 \ 2n+1 

n -2d-7/( Q) | Q ^n - S E J-=J_ // (Q )| Q ^n = o 

(65) 

In particular one sees that for quadratic interactions (local or nonlocal) the 
quantum evolution is formally identical to the classical one. In the configu- 
rational quantum cat we have a two-degree of freedom Hamiltonian 

oo 

i?(x, p, t) = # (x, p) + ff*(x, p) ^ ~ n ) > ( 66 ) 

n=— oo 

where both ff and are quadratic forms in the position and momentum 
operators. To write the functions explicitly, we define a 4-vector 



Q=m (67) 



and symmetric 4x4- matrices Bq and B k . 

Then the Hamiltonians Hq and if& are taken in the form 

H = ^QB Q , H k = ^QfifeQ (68) 
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The system with the Hamiltonian ( |66"D has four linear integrals of motion E5| 

I(t)=A(t)Q, (69) 
where the symplectic matrix satisfies the equation 

k(t) = A£5(t) , (70) 
with a 4x4-matrix E with identity 2x2-blocks 



1 
-1 



and 



B{t) =B + B k J2 S (t ~ n ) ■ 



The initial condition for the 4x4-matrix A(t) is 

A(0) = 1 

The Floquet solution to ( |7D| ) has the form 

A(l+) = e EB ° e SBfc 



A(n+) = A»(l+) = (■ 



and for n kicks 
The equation for the propagator is 
n +z 



)' 



(71) 
(72) 

(73) 

(74) 
(75) 



— i 



ij (x«,p«)+^(x( i ),p«) x; 

n=— oo 
oo 



n 



(76) 
11 = 0. (77) 



the vector-operators being 

-i -l 



,(i) 







P 



(i) 



i d 
2* l dX x 



did 



dm 2 8X± 







dXo 



-i 



" +^|(78) 



2 ^ax 2 



9 19 i d 

dX^ &v~2 M2 dx 



" 9 



-1 



<9z/ 2 
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and 



r( 2 ) 



■(!)* 



P 



(2) 



P 



(1)* 



The equation for the Fourier component G(k, /2, z7, t) used to compute the 
Lyapunov exponent is the same as Eq.fl76f) with 



d d 



dXx dX 2 
The equation (|76|) can be integrated by 

I! (X, /2, v, n) = II (X, /2 A , z?a, 0) 



(79) 



where the parameters /2a, z^a are expressed in terms of /2 and z7 by the matrix 
product rule 

(?a,/Za) = (^yu) A _1 (w+), (80) 

where A _1 (n + ) is given by ([75]). 

For the Hamiltonians Hi and of Eqs. (|58D and (|59| ) the matrices So 
and .Bfc are 



B n 



for i^x and 



5 n 



/ 1 













1 




















V 1 










/ i 








o\ 





1 




























0/ 



Bi. 



B, 











o\ 








1 








1 








\o 








0/ 


/o 








1 \ 








1 


1 





1 








V 1 


1 





0/ 



for H 2 . 

For the model [p]] where only the kick contributions are kept in the Hamil- 
tonian one has 

B = 



and 



Bt 



ln(l + u) 

uj + 2 



C(u) 
£M 
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the 2x2- matrix being 

-U) 



2(l+aj 

2uj on 



with oj = (1 + 

For this model the matrix A~ 1 (n + ) in (ISO) is given by 




where 

f 



For all three models the equations ([76]) are the same for classical and quantum 



motion. Therefore the Lyapunov exponent must be the same in the classical 
and quantum cases. In particular, as is known from the classical case0], 
there is a positive Lyapunov exponent, namely 

A = \nu 2 



5.4 The standard map 

The standard map is a case where the phenomena of wave function localiza- 
tion is believed to have a taming effect on chaos. The Lyapunov exponent 
analysis gives a characterization of how this taming effect comes about. 
The Hamiltonian is 

2 00 

H = P - + 1C0 s{q) Ht-nr) (81) 

n=— oo 

the configuration space being now the circle, q G S . This system describes a 
particle rotating in a ring and subjected to periodic kicks. It has been exten- 
sively used in studies of quantum chaos ( [fKJ - |)0J ) and has even been tested 



experimentally with ultra-cold atoms trapped in a magneto-optic trap |pl | . 
From (pH|) the equation to be solved now is 



dG 
~dt 



dG °° ( H \ 

S Ht-nr)sm(-uj{G{l,fi + l,u,t)-G{l,fi-l,u,t)} 

n=—oo ^ * 

(82) 
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where we have specialized to the value k = 1 because this is the only k— value 
needed to compute the Lyapunov exponent (fi3|). Notice that we have used 
here the same tomographic transformations that were described in Section 
2 for functions on the line. This is justified by considering all functions as 
defined not in S 1 but in the suspension of S 1 . 

From (jS2| ) one sees that between any two kicks the function propagates 
freely, namely 

G (1, /i, v, t ) -> G (1, fJ,, i/, ti_) = G (1, /J,, v + fj,r, t ) (83) 

and at the time of the kick a quantity is added that is proportional to a finite 
difference (in //). 



G (1, ix, v, t l+ ) = G (l, /i, !/, h_)+lf (y) {G (1, /i + 1, !/, *!_) - G (1, m - 1, *i 



7 

2 

(84) 

where, for the classical case 

f(u) = u (85) 

and for the quantum case 

/(") = ! Bin (86) 

To compute the Lyapunov exponent we need the evolution of the deriva- 
tives (in n and v) of G at ii — v — 0. From ( |83"D and ( |8^ ) one obtains the 
following iteration for the derivatives 

G^ (1,0,0, t+ 1) = G( 2 )(l,0,0,t) + G( 3 )(l,0,0,t) 
G( 3 )(l,0,0,t+1) = G( 2 )(l,0,0,t) + |(G(l,l,r,t)-G'(l,-l,-r,t)) 

(87) 

Let us consider first the classical case (h — , / (u) — v and 7 > 0). Let 
also r = 1 and go = Po — in the initial condition (pEBp. Then, one obtains 
the following recursion for the derivatives of G at fi = v = 0. 

G( 2 ) (1,0,0, n+ 1) = G( 2 ) (1,0,0, n) + (1,0,0, n) 

(1,0,0, n+ 1) = 7 G( 2 )(l,0,0,n) + (l + 7 )G( 3 ) (l,0,0,n) 

which has the solution 

G( 2 )(l,0,0,n) = A n {z)v x +B n {z)v 2 

G( 3 )(l,0,0,n) = C„(z)^i + J D n (z) U2 1 ' 
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with z = 2 + 7 and 

A n (z) = f4-i(f)-f4- 2 (|) 

C n (z) = C4(f)-2[/ n _ 1 (|) + C/„_ 2 (f) 

D n (z) = U n (l)-U n ^) 



(90) 



sin((n+l) < 
sin(cos" 

For the Lyapunov exponent one obtains in this case 



where U n ( Chebyshev polynomial. 



A = In 



(91) 



a result similar to the harmonic kicks on the line. One sees that as long 
as 7 > the exponent A in Eq.([H]) is always positive. This results from 
the choice made for the phase space point (po = qo = 0) where the marginal 
distribution receives the singular perturbation (|34j). If instead we had chosen 
(po = and go = ^) hi the initial condition ([45]), one sees easily by a change 
of coordinates in the Hamiltonian that this is equivalent to replace 7 by 
—7. Then the Lyapunov exponent A in Eq.(^l|) is positive only for 7 > 4. 
As discussed at length in the next section, this only means that it is the 
phase space point (po, go) that defines the measure for which the Lyapunov 
exponent is computed. Hence, for the measure that supports the hyperbolic 
point (po = , go = 0) the exponent is always positive, whereas for sufficiently 
small 7 > the exponent for the measure that supports the elliptic point 
(p Q = , g = 7r) is negative. 

For the quantum case (H 7^ 0) let us consider an initial condition G (1, /i, u, 0) 
fi + v (corresponding to po — , go = 0, V\ — i>2 = 1) and r = 1. Accord- 
ing to Eq . (jS7D , all one needs to compute the Lyapunov exponent is the time 
evolution of G (1, 1, l,t). For this purpose we set up a matrix recursion for 
the evolution equations (|83"1 - |34"D . Define the following matrices 

\ / 1 1 \ / 1 1 \ 

M = I 1 ; M+ = 1 ; M_ = 1 (92) 
) \) \ \ 1 I J \-l -1 -I J 

a \ 

and vectors ( (3 where a counts the number of /x's, (3 the number of v's 
7 / 

and 7 is a simple number. Then with Tr denoting the sum of the elements 



23 











x = 


1 


|;«b=| 






V 1 







in a vector and 

/ 1 \ / \ 

(93) 

the initial condition is G(l,l,r, 0) = Tr(x ) and the function / (z/) = 
f(Tr(y )). 

On arbitrary functions of 3- dimensional vectors, the operators Kq, K + , K_ 
act on the arguments by the matrices Mo, M+, M_ 

K i g(x)=g(M i x) (94) 

Then 

G (1, 1, 1, n)=Tr[ (k q + |/(y ) [#+ - q)"x } (95) 

where it is understood that the power of the operator is fully expanded before 
the Tr operation is applied to each one of the vector arguments. When this 
expansion is made, one obtains an expression of the form 



(n , k<(n — k)'. 

k=l i=l y 

(96) 

The products of M matrices in the arguments of / (•) contain a variable num- 
ber of factors, from 1 to k. However for each term, a different combination 
of products will appear. For h ^ the function / (v) is proportional to a 
sine and, if ^ is irrational, the coefficient of each •y k behaves like a sum of 
random variables of zero mean. Therefore each coefficient averages to zero 
and 

£(1,1,1,71)^71 + 2 (97) 

Large fluctuations are however to be expected in view of the large number 
of terms in the sums for large n. From (|5TD, the result ( |9"7| ) now implies 

G( 2 )(l,0,0,n) ~ 1 + n (1 - 7 ) + \ + \ 

G( 3 )(l,0,0,n) ~ l + 7 ( (ra+1) 2 ( " +2) - 1 ) 

For large n, logG^ 2 '' (1,0, 0,n) ~ 31ogn and XogG^ (1,0, 0, n) ~ 21ogn and 
the Lyapunov exponent vanishes. 
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The situation we have been studying (po = Qo = in the initial pertur- 
bation) corresponds to the (hyperbolic) case where the classical Lyapunov 
exponent is positive for any 7. We see here clearly the taming effect of quan- 
tum mechanics on classical chaos and its dynamical origin. It results from 
the replacement in the evolution equation of the linear function / {y) = u by 
/ (v) = I sin (f^). This in turn is a consequence of the replacement of the 
classical Boltzman equation by the quantum evolution equation fl2"2|) , or in al- 
gebraic terms, by the replacement of the ordinary product by the Moyal-Vey 
product in the non-commutative quantum phase space. 

In this model, the origin of the taming effect of quantum mechanics on 
classical chaos, is traced back to the existence, in the h— deformed equation 
flUD, of infinitely many terms in the series which add up to a bounded func- 
tion in v. How general this mechanism is, for other quantum systems, is an 
open question. In any case the taming effect of quantum mechanics, obtained 
here for the standard map, is more accurate than previous discussions of the 
same system, because it refers to the behavior of the Lyapunov exponent 
rather than to indirect chaos symptoms, like the energy growth or diffusion 
behavior. 

6 Remarks and conclusions 

1- The method developed in this paper, for the quantum Lyapunov expo- 
nents, provides a fairly unambiguous construction of these quantities, in the 
sense that classical and quantum exponents have the same functional form. 
The difference lies only on the time-evolution laws for the propagators. 

The dynamical evolution laws of the marginal distributions, obtained 
by the tomographic map, are apparently more complex that the familiar 
Schrodinger equation. However, for the computation of the Lyapunov expo- 
nents, they provide a fairly efficient computational scheme. 

Quantum mechanics is widely believed to have a taming effect on classical 
chaos. However, most discussions are of a qualitatively nature and fail to 
identify the conditions under which the taming effect is expected to occur 
and those in which it will not occur. This is a very relevant question in view 
of the fact that for local quadratic potentials, the quantum behavior differs 
very little from the classical one and genuine examples of quantum chaos 
with bounded configuration space are known, like the four-dimensional or 
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configurational quantum cat [0] j£| . 

In the standard map, studied in Sect. 5, it is clear that the suppression of 
chaos is directly related both to the nature of the potential and the analytical 
structure of the series in Eq. fl44]) . This operational series, when acting on the 
potential, convert an unbounded function into a bounded function in v (the 
symplectic parameter conjugate to p) . The structure of the series corresponds 
to the structure of the Moyal bracket and the way the quadratic potential 
(and presumably other polynomial potentials) avoid the suppression effect, is 
by truncating the action of the Moyal bracket to a finite number of cocycles. 

The fact that, for non-polynomial interactions, all derivatives of the po- 
tential intervene in the quantum evolution, means, by an analyticity argu- 
ment, that the future evolution of any local perturbation depends strongly 
on what is going on at all other points. This interference between quantum 
"trajectories" is probably the decisive factor that determines the nature of 
the quantum modifications of classical chaos. 

2 - A second important question concerns the support properties of the 
Lyapunov exponents that have been constructed. In classical mechanics, 
Lyapunov exponents are ergodic invariants. That means that they are defined 
in the support of some measure. In the construction (both classical and 
quantum) developed in Sect. 3, the Lyapunov is obtained from a singular 
perturbation of the X— coordinate at the point fiqo + upo for each pair (/z, v). 
For the classical case, the interpretation is clear. From the point of view 
of measures in phase-space, it means that one is constructing the Lyapunov 
exponent that corresponds to the measure whose support contains the point 
(?o,Po)- 

Measures on classical phase-space may be interpreted as measures on 
the joint spectrum of the (commuting) operators q and p. Therefore in the 
classical case a measure \i plays the double role of a probability measure 
in phase-space and a spectral measure for the dynamical operators. For 
the quantum case, however, q and p do not commute and there is no joint 
spectrum for these operators. Then, instead of one measure playing a double 
role we have two: 

- One is the state that is perturbed. This is the analog of the classical 
probability measure, because states are the non-commutative analogs of Borel 
measures. 

- The other is the spectral measure of the operator X, the perturbation 
acting, for each pair (qo,po), at the point p,qo + vpo of the spectrum. 
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In conclusion: the interpretation of the quantum Lyapunov exponent as 
an ergodic invariant requires two measures, a state and a spectral measure. 
In the classical case the two measures coincide. 
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